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ABSTRACT 

We compute the combined two and three loop order correction to the 
spin-spin correlation functions for the 2D Ising and q-states Potts model with 
random bonds at the critical point. The procedure employed is the renormal- 
isation group approach for the perturbation series around the conformal field 
theories representing the pure models. We obtain corrections for the corre- 
lations functions which produce crossover in the amplitude but don't change 
the critical exponent in the case of the Ising model and which produce a shift 
in the critical exponent, due to randomness, in the case of the Potts model. 
Comparison with numerical data is discussed briefly. 
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In severals works, the critical properties of the two dimensional Ising and Potts models 
with quenched random bonds has been investigated in a numerical and theoretical way. 
For the Ising model, it seems that the effect of impurities changes the divergent form of 
the specific heat near the critical point C(t) ~ Inln while the critical exponent 
of the < oo > correlation function is unchanged [fj, |3|. The renormalisation group (R.G) 
equations for the coupling constant as well as the correction for the < ee > correlation 
function has been established to the second loop order [[|. Recent numerical data gives 
interesting behavior in the crossover region for the < oo > correlation function || . These 
results suggest a deviation of such a correlation function in comparison to the pure case, 
deviation which is not predicted in previous perturbative computation up to the first 
order . We will extend in this letter the computation of this correlation function up to 
the third order. We found that there is a multiplicative correction in < oo > (compared 
to the pure 2d Ising model). We then compare these corrections with the ones of ||. We 
also consider the random Potts model, for which we introduce a perturbed central charge 
(c = n+ e )j which can be seen as a short distance regulator^]. We obtain the correction for 
the critical exponent of the < oo > up to the third order. The choice of the conformal 
field theory which should represent exactly (in a non perturbative way) the Potts model 
with random bonds at the critical point remains open and a complete non perturbative 
characterization of the critical exponents is needed. 

Our starting point is the Hamiltonian at the critical temperature : 

So = £ (-(3J )S Si , Sj (1) 

<i,j> 

where Jo is the coupling between bonds and s, are the Ising or Potts spin variables. The 
unperturbed partition function is obtained by taking a sum over the spins configurations : 

2o = E«" So (2) 

Si 

Near the critical point, the addition of a position dependent random coupling constant 
gives the following effective Hamiltonian expressed in a continues way : 

S = So + J d 2 z m(z)e(z) 
x Such a procedure was already used by A. Ludwig in the computation of the thermal exponent 0. 
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Here e(z) is the local energy operator of the critical pure model. Using the replica method 
f|], one can find the free energy by averaging over a distribution of the random function 

m(z) : 

Z N -1 



log(Z) = lim 



TV^O N 

with : 

TV . TV 

- S 0,a ~ / d 2 Z m(z) ^ £ a{z) 

Z N = J2e 0=1 a=1 (3) 

Si 

The average of Z N with an arbitrary distribution for m(z) gives an effective Hamiltonian 
containing all the cumulants of the probability measure. However, it can be shown by 
power counting that for the q-state Potts model all but the two firsts cumulants are 
irrelevant for g< 3 @j. The effective Hamiltonian can then be obtained by averaging the 
partition function with a Gaussian distribution : 

W= I l[dm(z)Z N e-^ {miz) - mo)2 

z 

which gives : 

TV , TV , N 

-J2 S o," +9o d 2 z £a{z)e b {z) - mo d 2 z Y £ a{z) 

~%n = e a=1 a ' b=1 a=1 (4) 

Si 

mo represents the deviation from the critical temperature of the model with random 

TV 

bonds and we are interested in the limit mo — > 0. Inside / d 2 z >J e a (z)eb{z), the terms 

a,b=l 

with coinciding replica number (a = b) can be omitted because they produce, under an 
operator product expansion (OPE), only divergent constant terms which can be removed 
by subtraction and irrelevant operators. Diagonal replica number terms will then be 
discarded in the following and the partition function can be written like 

TV . TV . TV 

-J2 S o,* + 9o / d 2 z J2 £ a(z)£b(z) - m / d 2 z Y £ a(z) 

^ e a=l J a^b J 0=1 (5) 

Si 

Correlation functions are calculated in a perturbative way around the conformal invariant 
action Sq. So a correlation function of some local operator O is expanded like : 

< O(0)O(R) >=< O(0)O(R) > + < SiO(0)O(R) > +- < S]O(0)O(R) > + • • • 



<>o means an expectation value taken with respect to Sq and 

. r N . N 

Sj = J d 2 z H^z) = 9o J d 2 z e a {z)e b (z) - m J d 2 z E £ *( z ) (6) 

a^b a=l 

The integrals of correlation functions involved in the calculation can be performed by 
analytic continuation with the Coulomb-gas representation of a conformal field theory 
H where the central charge is c = | + e. The e term corresponds to a short distance 
regulator for the integrals. In addition, we also used an I.R. cut-off r. The result is then 
expressed in a e series with coefficients depending on r. Then the limit e — * corresponds 
to the Ising model while the Potts model is obtained for some finite value of e. We recall 
here some notations of the Coulomb-gas representation for the vertex operators . The 
central charge c is characterized by a 2 + = = f + e with 

c = 1 - 24al ; a ± = a ± \j al + 1 (7) 

a + a_ = —1 

Note that for the pure 2D Ising model a\ = | and c = |. The vertex operators are 

V nm {x) = e ia ™*W ( 8 ) 
with 6( free scalar field and a nm defined by 

a nm = - (1 - n)a^ + ]- (1 - m)a + (9) 

The conformal dimension of an operator V nm (x) is A nm = —anma nm with 

1 1 

cw = 2a - a nm = - (1 + + - (1 + m)a + (10) 

The spin field can be represented by the vertex operator V PtP _i and V\p for the energy 
operator e. Note that in the Ising case the a operator could also be represented by the 
Vix operator (since both operators coincide in the limit e — > 0) and the final result turn 
out to be independent on which representation we are taking [0. 

We now turn to the computation of the renormalisation of o produced by the operator 
product 

k k N N N 

"° n J d 2 Xl d 2 z E £ ai (^K(^)E^) -+9aX(r,e) J d 2 zY.°a{z) (11) 



k\ 



The asymptotic behavior of the correlation function can be calculated with the help of 
the R.G. equations. A straightforward calculation yields a null contribution at the first 
order. The next two orders are 



_2 N N N 



2 

and 



J d 2 xd 2 yd 2 zJ2e a (x)e b (x)J2e c (y)e d (y)J2a e (z) (12) 



\x-z\,\y-z\<r a ^ b c ^ d 6=1 



3! 



d 2 xd 2 yd 2 wd 2 z Y e a (x)e b (x) Y £c{y)e d {y) Y e e (w)e f (w) Y & 

\ x - z \,\y-z\,\w-z\<r a + h c + d e ^f 9 



9\ Z ) 



(13) 

The renormalisation constants appearing at each order can be calculated by projecting 
(|l2l) and (13) over a (oo) giving respectively the following integrals : 

2g%(N - 1) J d 2 xd 2 y < a(z)e(x)e(y)a(oa) >< s(x)s(y) > (14) 

\x— z\,\y— z\<r 

4(N - 1)(N - 2)gl J d 2 xd 2 yd 2 w < a(z)e{x)e{y)a(oo) > x 

\x—z\,\y—z\,\w—z\<r 

< e(x)e(w) >< s(y)e(w) > (15) 



Using now the Coulomb-Gas formalism, the second order integral (H) can be rewritten 
like : 



2 g 2 (N - 1)N 2 J d 2 xd 2 y J d 2 u < V^ T (0)V 12 (x)V 12 ( 2 /)V p , p _ 1 (oo) J + ( M ) > \x - y\~ 2A ° 

\x\,\y\<r 

= 2g 2 (N-l)N 2 J d 2 xd 2 y J d 2 u \x\ 2a \y\ 2a \u\ 2b \x - u\ 2c \y - u\ 2c \x - y\~ 2 ^ (16) 

\x\,\y\<r 

with a = 2(2ao — a PtP -i)a\ 2 , b = 2(2ao — a PtP -i)a + ,c = 2a\ 2 a + and d = 2a.yio.vi- Here 
J+(u) is the screening charge operator 

J + ( u ) = e ia+ ^ u) (17) 

and N 2 is a normalization coefficient involved in replacing the correlation function < aeeo > 
by its Coulomb-Gas representation. At the third order, we have : 
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4(N-l)(N-2)g 3 N 2 x 
d 2 yd 2 zd 2 w I d 2 u < V^{Q)V 12 {y)V X2 {z)V p , p -x{oo)J + {u) > \y - w\- 2 ^\z - w\~ 2 ^ 



\x\,\y\,\w\<r 

= A(N-l)(N-2)g 3 N 2 x (18; 
J d 2 yd 2 zd 2 w J d 2 u\y\ 2a \z\ 2a \u\ 2b \y - u\ 2c \z - u\ 2c \y - z\ 2d \y - w\~ 2 ^\z - w\~ 2 ^ 

\x\,\y\,\w\<r 

The details of the calculation of these integrals (and for the normalization constants N 2 ) 
will be presented elsewhere 0. For the calculation of N 2 , see also |§. The final results 
for each of these integrals are 



-3(JV - l)<? Ve 



6e 



r 2 (-|)r 2 ^' 



1 a. 9 - 2- ^— 

n(-i)P(-i; 



for (HI) and 



3^3 



-9e~ 



12(JV-1)(JV-2)^tt 



9e 



1 + 



4 y 2 {-\)Y 2<v 



3r 2 (4)P(-| yj 



for (Jf5|). The a field gets then renormalised as 

a -> a(l + A 2 # 2 + A 3 #j + ...) = Z^tr 



and the constant Z CT satisfies the following equation 



dln{Z a ){r) 



-3(N -l)g 2 n 2 er- 6t 



dln(r) 

+12(N- l)(N -2)g*ir 3 r- 9t 



r 2 (-f)r 2 ^ 



1 + 2 



1 + 



r 2 (-|)r 2 H) 

4 r 2 (-|)r 2 (|) 
3r 2 (-|)r 2 (-i; 



(19) 



A' 



In the same way, the operator e a (x)£b(x) gets renormalised by the operator product 
of e operators coming from the perturbation term. This gives us the relation between 



the bare and the renormalised coupling constant : 



,3c 



g{r) = r *g + Att(N - 2)—g 2 + O(g 3 ) 



which produce the following R.G. equation 

dg 



dln(r) 



3eg + 4:7r(N - 2)g 2 



(20) 



(21) 
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Replacing g with the renormalised coupling constant g obtained from eq.([20D, the equa- 
tion for Z„ becomes : 



dln(Z a {r)) 



-3{N- l)g 2 (r)nh 



1 + 2 



r 2 (-pr 2 (|) 

r2 (-D r2 (-D 



dln{r) 

+ A(N -1)(N -2)ttVW ( 22 ) 
which is valid up to the order g 3 , g 2 e. 

The Callan-Symanzik equation gives the form of the correlation functions 

9« 



< a(0)a(sR) > a , so = e <* s ' 2 ^ < a(Q)a(R) > r , ff(s) (23) 

where we used the notation : = 1h{g) ', g(a) = 9o and g(s) is defined by 

/ (3(g)dg = ln(s); r = sa, and a is a lattice cut-off scale. In eq.(|23|), i? is an arbi- 

90 

trary scale which can be fixed to one lattice spacing a of a true statistical model. The 
dependence in s of < a(0)a(a) > r ,g{s) wm then be negligible assuming there isn't inter- 
actions between distances smaller than a and so it reduces to a constant, s will then 
measure the number of lattice spacing between two spins in < a(0)a(sR) >. In the 
following, we adopt the choice a — 1. 

The quenched model will be obtained in the limit iV — > 0. In the case of the Ising 
model, e — > and eq.(pl|) (with e = N = 0) gives the renormalised coupling constant 

^ ^ 1 + 8ngoln(s) ^ 
We can also compute the exponent in (|23|) 

?*-7'S*- : r(«M , -«8> w 

SO ^ Vay 90 

from which we can deduce the asymptotic form of the < oo > correlation function : 

Since we are considering the case where go is small, g(s) will also be small, and the 
exponent in (p6|) can be expanded. Using eq.fl2~3|) we obtain the final result : 

< >-~ f 1 + ^° f 1 " ( i +8 ,U(.) )')) (27) 



We can see that the third order correction increase the correlation function in the cross- 
over region (ln(s) ~ while at large s, as expected, the exponent is unchanged. 

We now turn to the case of the 3-state Potts model where a 2 , = f , e = — ^. Here 

t 5 ' 15 



lh{9c) + and eq.(0) will be (with N = 0) 

(3(g) = -3eg - Sirg 2 (28) 

Contrary to the Ising model case (e = 0) where /3(g) had an I.R. point at g — 0, here, 

because e is negative, g = become an U.V. fixed point and the new I.R. critical point 

9{s) 

is located at g c = — J| + 0(e ). So, the integral J j$N<ig is dominated by the region 

go (3> 



where g ^ g c and J ^jj^-dg « 7/ 1 (g c )/n(s). We thus obtain 



.</» 



< a(0)<r(s) > ff0 ~ s -( 2A - 2 7a gc )) (29) 



which show us that the new exponent for the dimension of the <r operator is 

r 2 (-Dr 2 ^ 



A' CT = A CT - 7fc (0 c ) = A CT - 37r^ 2 e 



Inserting the expression of g c , this becomes 
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3^ V6' 



r 2 (-|)r 2 (-|; 



A« (so) 



27 p^-^r 2 ^ 

a ^ a --1 p(-|)p(4) - 3 + °<- 4 » < 31 > 

The critical exponent of the spin-spin correlation function of the pure model is 2 A^ = ^ 
while the new exponent produced by the randomness is 2A^ = ^ + 0, 00264 = 0, 26931. 
We see that unlike for the case of the thermal exponent [f|], the first correction appears 
at the third order. The change of the exponent is of 1% of the total magnitude. 

Recently the numerical simulations of the random Ising model which measure directly 
the deviation of < oo > from the pure Ising model at the critical point has been per- 
formed ||. These measurements were made for disorder such that 87rgo ~ 0.3 |f|, [10 



Deviations predicted by our computations are very small. They correspond to 0.1%. The 
deviations obtained in numerical simulations are around ten times larger, and they are 
of opposite sign, i.e. correspond to a decrease of the spin-spin function with distance 
r. In |J, it has been checked that this decrease corresponds, within the accuracy of 

8 



the measurements, to a factor function of the ratio r/L, F(r/L), r being the distance 
between the spins and L is the lattice size. So they correspond to finite size effects, 
being different for perfect and random models. We would suggest, on the bases of our 
calculation of the purely r dependence of the spin-spin function on an infinite lattice, 
that numerical deviations will continue to be plotted by the same curve F(r/L), if one 
measures < oo > for different lattice sizes as it has been done in ||, until the accuracy 
reaches the value of the r-deviation which we calculated here. Only then the curves for 
different L would split. Probably it would be much easier to observe the deviation in the 



case of the random 3-states Potts model. There, according to our calculation eq.(^T|), it 
is the magnetization exponent which would decrease by about 1%. 
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